ROMEl /1446-06 
IFIC/07-22 
hep-ph/ yymmnnn 



The Two Loop Crossed Ladder Vertex Diagram 
with Two Massive Exchanges 

U. Agliettfl 

Dipartimento di Fisica, Universitd di Roma "La Sapienza" and INFN, Sezione di Roma, 

1-00185 Rome, Raly 

R. Boncianill, 

Departament de Fisica Tedrica, IFIC, CSIC - Universitat de Valencia, 

E-46071 Valencia, Spain 

L. Grassil, 

Dipartimento di Fisica dell' Universitd di Bologna, and INFN, Sezione di Bologna, 

1-40126 Bologna, Raly 

E. RemiddH 

Dipartimento di Fisica dell' Universitd di Bologna, and INFN, Sezione di Bologna, 

1-40126 Bologna, Italy 



Abstract 

We compute the (three) master integrals for the crossed ladder diagram with two 
exchanged quanta of equal mass. The differential equations obeyed by the master in- 
tegrals are used to generate power series expansions centered around all the singular 
(plus some regular) points, which are then matched numerically with high accuracy. 
The expansions allow a fast and precise numerical calculation of the three master 
integrals (better than 15 digits with less than 30 terms in the whole real axis). A con- 
spicuous relation with the equal-mass sunrise in two dimensions is found. Comparison 
with a previous large momentum expansion is made finding complete agreement. 
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1 Introduction 



Electron-positron linear colliders of next generation with a large energy, E = 0{lTeV), 
and a large luminosity, L = O (10^^cm~^s~^), will produce an incredible amount of accurate 
data [IJ. Theoretical work in the same direction, i.e. pointing to accurate computations 
of electroweak cross sections and decay widths, is therefore mandatory. That means, in 
practice, to go beyond the standard one-loop approximation and address the control of the 
two-loop quantum corrections. While in QCD multi-loop computations are often feasible 
by setting to zero the small quark masses simplifying then the problem to a massless one, 
in the electroweak case setting to zero the masses of the W, Z, Higgs or top quark is not 
possible. One has to face the general problem of computing massive multiscale Feynman 
diagrams. In many cases, the full computation can be drastically simplified by finding 
a small expansion parameter t], such as the ratio of particle masses and/or kinematical 
invariants: 

. = (1) 

The heavy algebra entering the higher orders in rj can be dealt with using standard alge- 
braic programs such as Form [2] or Mathematica [3J. In problems involving many different 
scales or, equivalently, many different dimensionless ratios, that is probably the only viable 
analytical way. Another strategy is that of numerical computation of the integrals "from 
the beginning" , such as for example the numerical computation of loop integrals written in 
terms of Feynman parameters [1] . This approach is not easy to implement when a subset of 
particles is massless, because of troublesome infrared effects. A different strategy assumes 
the dominance of the logarithmic terms over constants and power-suppressed terms in the 
high-energy cross-sections, i.e. of terms of the form [5] 



where g is the weak coupling and m is the W or mass. The logarithmic terms are 
computed directly by means of asymptotic expansions P, [7]. This strategy is similar in 
spirit to the one used in QCD to construct for example shower Montecarlo programs [8]. 
The logarithmic terms are certainly the dominant ones in the formal limit s oo (i.e. 
s S> m), but at a fixed energy non logarithmic terms may in general have a significant 
numerical effect. In cases in which one or (at most) two dimensionless quantities of order 
one are involved, an exact analytic computation may be useful. A model computation for 
higher-order electroweak effects is the two- loop form factor in the degenerate mass limit. 



In [9], [To] two of us have computed the master integrals of the two-loop electroweak form 
factor for all diagrams with the exception of the crossed- ladder with the exchange of two 
equal mass quanta — such as a pair of W^s and/or Z's (see fig. 1). These master integrals 
are naturally expressed in the basis of the harmonic polilogarithms (HPLs), transcendental 
functions verifying algebra properties which largely simplify the transformations [TT] . Some 
master integrals involving two or three massive propagators have required a generalization 
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of the harmonic pohlogarithms basis, involving basic functions with square roots — the 
so-called generalized harmonic pohlogarithms (GHPLs) [lOj. 

In this paper we compute the last missing piece in order to obtain an exact analytic 
expression of the two-loop electro-weak form factor: the crossed ladder diagram with two 
equal-mass quanta exchanged. By standard reduction techniques, we find that this topology 
involves three master integrals. At variance with respect to the planar topologies and the 
non-planar topology with at most one massive exchange, it is not possible to represent 
any of these master integrals in terms of harmonic polylogarithms or generalized harmonic 
polylogarithms: elliptic integrals are involved, i.e. integrals of square roots of polynomials 
of degree three or four [12]. For that reason, we use a different method, already applied 
in [131 [E] to evaluate the equal-mass sunrise diagram. The differential equations obeyed 
by the master integrals are used to generate power-series expansions centered around all 
the singular points of the differential equations themselves. It may also be convenient 
to generate auxiliary expansions around regular points, which are in principle completely 
arbitrary. The series expansions centered around two different points are then matched 
numerically in a point belonging to the intersection of the respective domains of convergence. 
In general, we are not able to find a closed analytic expression for the master integrals, 
but only truncated series representations. The latter allow, however, a fast and accurate 
numerical evaluation of the master integrals themselves (better than 15 digits with less 
than 30 terms, typically, for the "accelerated" versions of the expansions). The logarithmic 
terms in the large momentum expansion of one of the master integrals (the basic scalar 
amplitude) have been obtained in [15j. Our expansion is in agreement with these results; 
we are also able to compute the power-suppressed corrections. 



2 Power Series Solution of the Differential Equations 
on the Mis 



In this section we outline the computation of the master integrals of the two-equal-mass- 
exchanged crossed-ladder diagram, which is described extensively in the next sections. As 
anticipated in the introduction, we compute the master integrals with a general semi- 
analytical method [12], which can be applied to arbitrarily complicated cases. The first 
steps (sec. H]) are common to the method used to compute the previous electroweak form 
factor diagrams [HI [ID]- By using standard reduction procedures, the problem is shifted 
to that of computing three independent scalar amplitudes, called master integrals (Mis), 
which obey a system of coupled first-order differential equations in the evolution variable 
X = —s/m^, where s = is the Mandelstam variable (sec. [5]). These equations are linear 
and not homogeneous. By a suitable choice of the master integral basis, we succeeded in 
triangularizing the system to a system of two coupled equations (let's say the first two 
equations on Fi and F2) plus a decoupled equation from the previous ones (let's say the 
third one on F3): 




(4) 



(5) 
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dx 



= A^i{x)Fi{x) + As2ix)F2{x) + A-isix)Fsix) + n^ix), (6) 



where the Aij{x)^s, the coefficients of the associated homogeneous system, are rational 
fractions in x, while the Qk{x)^s are known functions (logarithms, polylogarithms etc.). We 
then derive from the system of the two coupled ffist-order equations a single second-order 
equation in one of the two master integrals involved, let's say Fi (sec. |6]): 

r/^Fi dF, 

-4 + A{x)^ + B{x)F,{x) = n{x). (7) 

The resulting equation belongs to the Fuchs class, i.e. it has regular singular points only, 
including also the point at infinitj0 [12]. There are four such singularities, all located on 
the real axis, at: 

a; = 0; a; = 8; x = —1; x = oo. (8) 

The differential equation ([7j) is used to generate power series expansions centered around all 
these singular points (sec. [7j). Auxiliary expansions around (arbitrary) regular points may 
also be of utility. Each series is determined by the differential equation up to two arbitrary 
coefficients, which must be fixed by some initial or boundary condition. A qualitative 
knowledge of the behaviour of the master integral Fi[x) around x = is sufficient to fix 
exactly (in fact analytically) the two arbitrary coefficients of the expansion around this 
point. The coefficients of the series centered around x = (small-momentum expansion) 
are therefore completely determined in analytical way. The arbitrary coefficients occurring 
in the serieses centered in the remaining points are found by imposing matching conditions 
in suitable points belonging to the intersections of the respective domains of convergence. 
Assume that one knows the expansions of -Fi(x) around two different points xi and X2 

oo 

Fii.^) = X^OnH^^ ~ 2;i)" for |x-xi| < i?i , (9) 

n=0 

oo 

Fi(x) = J]a(f)(x-X2)" for |x-X2|<i?2, (10) 



n=0 



with |xi — X2I < -Ri + -R2, so that that the two circles overlap, and that the coefficients On ■* 
of the ffist series are completely known, while those of the second series are determined as 
(linear) functions of the ffist two by the underlying differential equation: 

a^^^ = hnttf'' + c„af^ + dn- (11) 



^ In general, by a singularity of the differential equation we mean a singularity in the coefficients A{x) 
or B{x). A point a;o is a regular singular point of a differential equation if its solutions can be written as 
singular factor ~ (x — xq)" or ~ (x — xq)" log(x — xq) multiplied by convergent serieses in a neighborhood 

of Xq. 
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By taking a point x in the intersection domain, the matching conditions give a hnear system 
on the unknowns a^^^ and a^^. 

oo oo 

n=0 n=0 

00 00 

^a^^^ (x - xi)" - (x - ; (12) 

n=0 n=0 

00 00 

Oq^^ n 6„ (x — X2)"^^ + af^ ^^nc„ (x — X2)"^^ = 

n=l n=l 

00 00 

77, a^^-* (x — xi)" — ndn{x — X2)"~^ . (13) 

n=l n=l 

Once aQ^"* and a^^"* are known, one determines the higher-order coefficients a'n' for n > 1. 
Instead of matching the derivatives of the serieses, one can also take a second matching 
point x' . 

Once the first master integral Fi has been determined, one plugs its expression into the 
first-order equation (jlj) of the 2x2 system involving dFi/dx and determines algebraically 
F2 (sec. [H])- One finally inserts the expression of the first two Mis Fi and F2 into the third 
differential equation ^ and determines F3 by quadrature. Let us remark that one can 
study all the analytical properties of the Mis (analytic continuation, behaviour close to 
thresholds, asymptotic properties, etc.) by looking at their power-series expansions. The 
latter also offer a suitable and powerful way for the precise and fast numerical evaluation 
of the Mis, for example with a Fortran routine. 

In sec. Owe discuss the relation between the differential equation on the master integral 
Fi for the crossed ladder and that one for the sunrise diagram with three equal masses 
in two space-time dimensions [13]. The "physical" origin of this connection is, at present, 
unclear. Finally, in section [10] we draw our conclusions. 



3 Threshold Structure 

The crossed ladder diagram with two equal-mass quanta exchanged has a rather simple 
threshold structure: there is a 2-particle cut on the internal massless lines for 

s > (14) 

and two 3-particle cuts on the internal massless lines and one of the massive lines for 

s > m^. (15) 

There is also a pseudothreshold for s < 4m^. The complexity of the differential equations 
and of the structure of the master integral cannot be guessed by looking at the threshold 
structure of the diagram. That is to be contrasted to the case of the equal-mass sunrise, 
where the threshold at s = 9m^ as well as the pseudothreshold in s = m? are identifiable 
as the basic sources of complexity. 
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Pi 



P2 

Figure 1: Feynman diagram for the annihilation of a pair of massless fermions with the 
exchange of two massive quanta with equal mass m. The thin lines represent the massless 
fermions, while the thick lines represent the massive quanta. The outgoing dashed line 
represents the probe (for instance a Z'). 

4 Reduction to Master Integrals 

By standard decomposition into invariant form factors and rotation of the scalar products, 
one can show that the computation of the two equal-mass crossed ladder diagram (see fig. [T]) 
is equivalent to the computation of the following independent scalar amphtudes: 

pm pn2 pna pn4 png pne '^^hV^k2, (16) 
1 -'2 -'3 -'4 -'5 6 

where D is the space-time dimension, the scalar product is defined as 

a ■ b = a ■ b — aobo, (17) 

the loop measure is 



r(3-D/2) 47r^/2' ^ ' 

with T{z) the Euler Gamma function. We consider a routing of the loop momenta and 
fcg which results in the following denominators: 

= kf + m^ (19) 

Pa = kl + m\ (20) 

^3 = {Pi - h)\ (21) 

Pa = (P2 - k^f , (22) 

= {Pi -h + k^f , (23) 

= {P2 + k,- k^f , (24) 

and the following irreducible numerator (scalar product): 

S = p2-ki. (25) 

The indices of the denominators are assumed to be all positiv^, > while the index of 
the scalar product can be positive or zero, r > 0. 



If < for some i we have a sub-topology in which line i is shrinked to a point. 
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The above independent scalar amplitudes, considered as functions of the integer indices 
rij's and r, are related to each other by integral identities obtained in general by means 
of integration by parts over the loop momenta, invariance under Lorentz transformations 
and symmetry relations coming from the particular mass distribution of the diagram under 
consideration [TB]. These identities can be used to express a given amplitude F as a linear 
combination of a set of reference amplitudes, with fixed indices, called master integrals (MI) 

N 

F (rii, ns, ng, n^, ng, Uq, r) = ^ q {ui, ns, ng, n^, rir,, rie, r) Fi, (26) 

i=l 

where the Cj's are known coefficients (rational fractions in x and D). The above reduction 
involves the solution of the integration-by-parts identities in some recursive way over the 
indices. Nowadays the most common algorithm is the "Laporta algorithm" [T7j, which 
we have used for the crossed ladder topology. That way we have been able to reduce the 
independent amplitudes to a linear combination of three master integrals: 

3 

F (ni, n2, ng, n^, n^, ng, s) = ^ q (ni, ^2, ng, n4, ng, ne, s) + • • • , (27) 

i=l 

where the dots denote contributions from master integrals with less than six denominators, 
i.e. from subtopologies, which are known [9l [TOl [18]. Let us choose the following basis of 
master integrals: 

F3(x;e) = [ V'^hV^k2. 

The first MI is the basic scalar amplitude, with linear denominators and no scalar products; 
the second MI contains the first denominator squared, while the third MI involves an 
irreducible numerator. By direct inspection, one finds that the above Mi's do not contain 
any pole in Z) — 4 because they are ultraviolet finite and, since the radiated quanta are 
massive, there is an effective infrared cut-off ~ m which renders the Mis infrared finite as 
well. Let us therefore set D = 4 from now on and omit this argument. 

5 The System of Differential Equations 

Once the reduction to master integrals has been achieved, the following step is the actual 
computation of the master integrals themselves. We use the differential equation method 
[T9] . which involves taking a derivative with respect to 

x = -^, (31) 
6 



(28) 
(29) 
(30) 



at fixed m^, of tlie master integral Fi. 



X 



m 

dx 



1 ^.m 

2 dpt 



1 ,m 



(32) 



q = Pi + P2 tlie momentum of the probe and = —s, s being the squared cm. energy. 
By taking the derivative inside the integral in eq. ( l32i) . various amplitudes are generated, 
which are reduced again to master integrals as discussed in the previous section. We then 
obtain in general a system of linear differential equations with variable coefficients, of the 
form: 

d ^ 

-Fi{x) = + Ni{x) {t = 1,2,3). (33) 



dx 



By changing the basis for the master integrals, the functions fij{x) and Ni{x) are trans- 
formed into new functions. A general basis for the Mi's involves a system of three coupled 
differential equations, which are equivalent to a single, third-order differential equation, 
whose solutions, as well known, are rather difficult to find. With the basis given in the 
previous section, the system comes out to be triangular with a 2 x 2 block. The coefficients 
of the associated homogeneous system read, in this basis: 



fuix) 



X 



/2l(x) = -— + 



+ 



1 



2x 3(x+l) 6(x-8) 



/3i(a;) = - — ; 

bx 



fuix) = -; 

X 



f32{x) 



X 



1 



X 



/i3(x)=0; 
f23{x)=0; 



fssix) = ; 

X 



(34) 
(35) 

(36) 



while the known terms, related to the subtopologies, are: 



Ni{x) 
N2{x) 



Nsix) 



0; 



16x 



^HiO,0;x)-^H{0,-l;x) + 



, ] Hir,0;x) 
^yx[A — X) 



IT 

T 



(37) 



A{x + 1) 
1 

"8(x-8) 
1 



-H(0, 0; x) - -H(0, -1; x) + , ^ 

9 ^ ' ' ^ 2 ^ ' ' ' 12^x{A-x) 



H{r,0;x) + — 



—H(0, 0; x) - -H(0, -l;x) , ^ 

36 ^ ' ' ^ 8 ^ ' ' ^ 3y/x{A-x) 



H{r,0;x) + — 



48x2 



^Hi-l, 0, 0; x) + H{r, r, 0; x) + x) 



(38) 
(39) 



The known terms are written in terms of generalized harmonic polylogarithms [lO]. Since 
fi3{x) = and /23(a;) = 0, F3 decouples from the system of the first two Mi's Fi and F2. 
It is therefore natural to compute Fi and F2 first. Once the latter are known, the third 
master integral F3 can be computed by quadrature — i.e. by integration of the first-order 
differential equation (l33l) with i = 3. 
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5.1 Master Integrals Close to Zero External Momentum 

Since pf = p\ = 0, 

q^ = 2pi-p2 (40) 

and the hmit 

X ^ (41) 

is equivalent to the limit pi ■ P2 —>■ 0. However, limit fHTj) also implies the stronger limits 
Pi —>■ and P2 ^ 0, which reduce the vertex topology to a vacuum one. That can be 
seen for example by introducing Feynman parameters Xi, X2 - ■ - x^ for the vertex amplitude 
and integrating over the loop momenta fcf and k!^. After that, the integrand can depend 
on relativistic invariants only, i.e. just on g^. One then performs analytic continuation to 
Euclidean space, where the limit ^ implies q'^ — > 0, i.e. P2 = —Pi- The diagram can 
then depend only on = and therefore must be equal to that one computed for = 0, 
i.e. to the corresponding vacuum amplitude. In particular, the Mi's reduce to the following 
vacuum amplitudes: 

FUx = 0) = [ ^V%V%; (42) 

J {kf + m^) [k^ + m^) kf kl [{ki - k2yf 

F2(x = 0) = [ i ^V^kiV^k2\ (43) 

^ ' J {kl + m^f{kl+m'')klkl[{ki-k2Yf 

Fs{x = 0) = 0. (44) 

The above integrals can be computed exactly by using the identities 

11 1 1 



kf kf + kf ni^ {kf + m^) 



(45) 



For our purposes, that is actually not necessary. By shrinking the massive lines in the 
infrared regions kf, k^ m?, we obtain: 

1 r^uv 1 

F,{x = 0) ^ — / — — V%V%; (46) 

J kf kl [{ki - k2f] 

1 r^uv 1 

F2{x = 0) ^ — / — — V%V%, (47) 

J klkl [{kl - k2f] 

where an ultraviolet cutoff Kuv ^ fnis assumed. The integrands above are invariant under 
the limit A — > in = 4, where A is introduced through the rescaling 

k^^ Xkf^, (48) 
fc^ ^ Afc^ (49) 

All that implies an over-all logarithmic soft divergenc^. There is also a logarithmic coUinear 
singularity for k^ cx k2 related to the denominator [{ki — ^2)^]^. As a consequenc^: 

Fi,2{x) ~ log^x, Fsix) ~ forx-^0. (50) 



There is not any soft sub-divergence, as can be seen by power- counting under the rescahngs 
A , ^ ^2 and k-^ ^ k-^ j kt^ ^ A k<2 ■ 

^ If one sets a; = from the very beginning, 1/(1? — 4) poles of infrared nature do appear in Fi and F2. 
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That implies in particular that no pole terms ~ l/x can appear. As we are going to show 
explicitly, this qualitative information is sufficient for the solution of the system at small 
external momenta. Finally, let us note that we can set m = 1 from now on without loosing 
any information (m is an over-all scale). 



6 Second-Order Differential Equation for Fi{x 

It is convenient to transform the system of differential equations for Fi and F2 into a single 
second-order equation for Fi only. The procedure is standard: one takes an additional 
derivative with respect to x on both sides of eq. (l33ll with i = 1; dFi/dx and dF2/dx are 
replaced by the r.h.s. of eqs. ( l33l) with i = 1 and i = 2 respectively; finally F2 is replaced 
by its expression in terms of dFi/dx and Fi, coming from the ffist of eqs. fl33l) . We obtain: 



Fi d Fi 

' (x) + A{x)—^{x) + B{x)Fi{x) + C{x) 



where: 



A{x) 
B{x) 

C{x) 



dx"^ 



3 1 

- + 7 

X X + 1 

^ + ^ - 



dx 



0, 



x^ 



+ 



5X 



1 



^ x-8^ 
4 

3(x + 1; 
7 



+ 



24(x 
1 



(51) 

(52) 
(53) 



1 



16x2 
9 

32x2 
1 

!2 ~ 



128x 18(x+l) 1152(x-8 
63 1 1 



256x 
7 



+ 



4(x+ 1) 
5 



+ 



256(x 
1 



8x2 32a; 24(x + l) 96(x - 8 



if(0,0;x) 

/J(0,-l;x) 
1 



X 



=i/(r,0;x) 



7 



+ 



5X^ 



64x 9(x+l) 576(x- 8) 



(54) 



Let us make a few remarks: 



1. the coefficient A{x) of -^{(x) contains simple poles, while the coefficient B{x) of -Fi(x) 
contains at most double poles. That is the necessary and sufficient condition on the 
differential equation to have regular singular points only 



2. the denominators entering the coefficients are: 



1 



X 



1 



1 



x + 1' 



4-x 



X — 



(55) 



The last denominator was not expected a priori^ as it does not correspond to any 
threshold or pseudothreshold of the diagram. The denominator 1/(4 — x) only ap- 
pears in the known function C(x) and it is related to the sub-topologies with the 
pseudothreshold in s = — 4m2; 
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3. the known term C{x) contains GHPL's (see appendix of weight at most two; 

4. the GHPL H{r,0; x), containing a single square-root basic function, has a coeffi- 
cient involving a square root 1 / a/ x{A — x) to ensure reality of the solution across the 
pseudothreshold located at x = 4. 

7 Power-Series Solution of the Second-Order Differ- 
ential Equation for Fi{x) 

The second-order differential equation (!5T|) for Fi is by far too complicated to be solved 
in a closed analytical form. We therefore look for solutions in the form of power series 
expansions centered around some points on the real axis. Let us recall that it is in any 
case necessary to consider the expansions around all the singular points, as the latter are 
by definition outside the convergence region of the ordinary expansions at regular points. 
Thanks to the differential equation satisfied by the considered functions, expansions around 
singular or regular points can be performed almost in the same way. 

Since the coefficients A and B vanish for a; = 0, —1, 8, we conclude that these points 
are singular points for the differential equation. Also x = oo is a singular point for the 
differential equation, as can be seen by changing variable to y = 1/x and studying the 
limit y — i> 0. 



7.1 Expansion Around x = — Small Momentum Expansion 

In this section we consider the expansion of the master integral Fi{x) around x = 0. Since 
the nearest singularity to the origin is located in x = — 1, we expect a radius of convergence 

Ro = 1, (56) 

i.e. convergence in the complex x-plane for 

|x| < 1. (57) 

For real x, that means: 

- 1 < X < 1. (58) 

Because of linearity, the general solution of the complete inhomogeneous equation can be 
written as the superposition of the general solution of the associated homogeneous equation 
— i.e. eq. with C(x) = 0, 

^^-^{x) + A{x) ^(x) + B{x) Fi^\x) = 0, (59) 
ax^ ax 

plus a particular solution of the complete equation: 

Fi(x) = Fi°\x) + Fi(x). (60) 
We will use systematically eq. (l60!l in the following expansions. 
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7.1.1 Homogeneous Equation 

Let us first consider the associated homogeneous equation. Since zero is a "singular regular 
point", we can look for solutions having the form of a singular function in x = 0, S{x), 
multiplied by a regular (i.e. convergent) power-serie^ : 

oo 

(x) = S{x) 5^A„a;^ (61) 

n=0 

where are numerical coefficients determined from the differential equation itself and 
from some initial or boundary conditions. The function S{x), giving the leading singularity 
for X — 0, is assumed to be of power-like form 

S{x) = x", (62) 

and solves the limit of the homogeneous equation for x ^ 0: 

S"ix) + -S'ix) + ^Six) = 0. (63) 

By inserting a solution of the form ( l62i) we obtain the indicial equation {a + 1)^ = with 
a double zero in a = — 1, implying two independent pre-factors of the form: 

SW = i!^. (64) 

X X 

The most general solution of the homogeneous differential equation is therefore of the form: 

oo oo 

Fi°\x) = J2 + logx ^ hnX"" , (65) 

n=—l n=— 1 

where we have absorbed a 1/x factor inside the series by defining: 

a„ = An+i. (66) 

By expanding the differential equation (!59l) around x = and substituting the series rep- 
resentation in eq. fl65l) . we can obtain recursively all the desired coefficients. The first few 
are: 



ao = 


1 3 

a_i o_i : 

4 8' 


bo = 




(67) 


Oi = 


5 33 , 

— a_i H o_i ; 

32 128 


bi = 




(68) 


02 = 


7 25 , 

a_i 0-1 : 

64 128 


62 = 


64 ' 


(69) 


as = 


173 2561 , 

a_i H o_i : 

2048 16384 ' 


^3 = 


173 , 

0-1 . 

2048 


(70) 



Let us make a few remarks: 

^ This case is to be contrasted to that of an "irregular singular point" , in which the singularity of the 
differential equation is so strong that no factorization of the singularity of the form is possible. In the 
latter case, one typically obtains series with an infinite number of negative powers or asymptotic (divergent) 
expansions. 
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• setting to zero for example a_i and all the higher-order coefficient vanish and we 
obtain the trivial solution: that is the homogeneity property; 

• there is a triangular structure: all the coefficients bi are proportional to the lowest- 
order one b-i, implying that by setting 6_i = we obtain a solution without the 
series with the logarithmic prefactor and with the simple pole in a; = only. The 
coefficients depend instead on both a_i and implying that by setting a_i = 
one obtain a solution containing both serieses; 

• A ffist independent solution can be obtained by taking for example a_i = 1 and 
6_i = 0, resulting in a function without logarithmic terms. A second independent 
solution can be obtained by taking a_i = and 6_i = 1, resulting in a function having 
both the power and the logarithmic terms. With a pictorial language, we may say 
that the logarithmic series "feeds" the standard one, while the vice-versa is not true. 

• The singularity of the differential equation at x = produces a logarithmic branch 
point in the solution at x = whenever 6_i 7^ 0. 



7.1.2 Complete Equation 

Let us now consider the complete equation fl5T]) . We have to expand around x 
the inhomogeneous term C{x), which is known (see appendix [B| : 



C{x) = ^ 



knX^ + log X ^ qnX^ + log^ X ^ r„x" . 



n=0 



n=0 



n=0 



The known term has a double pole in x = 0, multiplied also by a single or a double lo; 
of X. The expected radius of convergence of the multiplying series is one: 



R 



1. 



The explicit expressions of the lowest-order coefficients read: 



k2 

h 



8 32 ' 

23 7 
+ 



-71 



288 256 
3931 57 



-TT 



28800 2048 
1789247 455 



+ 



11289600 16384 



-IT 



?3 



19 

192 ' 
_ 671 _ 

^7680 ' 
38791 

430080 



r3 



1 

^32 ' 
7 

256 ' 
57 

^2048 
455 

16384 ■ 



also 

(71) 
;arithm 

(72) 

(73) 
(74) 
(75) 
(76) 



Let us look for a particular solution of flHT]) of the form: 

00 00 00 

Fii^) = Pnx'' + log X ^ g„x'' + log^ X ^ c„x" 



(77) 



n=-l 



n=-l 



n=0 
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By inserting the above form of the solution into the differential equation expanded around 
X = with the known term C{x) given by the series expansion in eq. (17T|) . we obtain for 
the coefficients: 

1 3 7r2 1 

= -4^-^-8^-^ + 32-16' 



Co 



Cl 



1 

32 



5 33 5 



= 32^-^+128^-^- 64 + 576' 



5 1 

- 32^-^ + 96' 
1 

64' 



7 25 13 2 1093 



64' 128' 1152 518400 
7 121 



C2 



64^ 17280 
13 



1152 



173 2561 57r2 649 



2048"^ 16384^ 576 635040 ' 
173 95 



93 = 77777779-1 

C3 -- 



2048^ 24192 
5 



576 



563 42631 407 ^ 3217 



P4 = -7777777^-1 - 777;;777779-l + 77777777 



8192' 327680' 57600 6720000 
563 1141 



94 = -^7^7779-1 

C4 



8192^ 432000 
407 



57600 

The coefficients Cj of the double-logarithmic terms are completely determined, while the 
remaining ones Pi and qi are fixed once two arbitrary coefficients, such as for example p_i 
and have been fixed. That is exactly the same arbitrariness that we have already 
found for the associated homogeneous equation. Since the general solution of the latter has 
already been found, we need only to find a particular, i.e. a single solution of the complete 
equation. Let us choose for example: 

p_i = and g_i = . (93) 
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With the following arbitrary choice, the particular solution does not contain any pole term 
and its numerical coefficients are completely fixed: 

Po = -^ + ^; go = -^; Co = l; (94) 

5 1 1 ^ ^ 

^^= 576 - 64' ^^ = 96' ^^ = "64' ^^^^ 

1093 137r2 _ _ 121 _ 13 

~ "518400 ^ 1152 ' 17280' ~ 1152 ' ^ ' 

649 57r2 _ 95 ^ _ 5 

~ 635040 ~ 576 ' ~ 24192 ' ~ ~576 ' ^ 

3217 40777^ _ _ _ 1141 _ 407 

~ "6720000 ^ 57600 ' ~ ~ 432000 ' ~ 57600 ' ^ ^ 

The general solution of the complete equation is therfore: 

oo oo oo oo oo 

^1 i^) = X] + log X ^ + ^ pnx"- + log X + log^ X ^ CnX"- . (99) 

n=— 1 n=— 1 n=0 n=0 n=0 

To uniquely determine the solution, one has to impose some boundary or initial condi- 
tions. As already discussed in the previous section, Fi{x) can have at most a logarithmic 
singularity for a; — > 0, implying that the coefficients of the power singularities must vanish: 

a_i = 0, 6_i = 0; (100) 

and the solution of the homogeneous equation to be selected, reduces to the trivial one. 
The particular solution Fi of the complete equation that we have chosen is therefore the 
expansion of the master integral Fi as defined by eq. fl28l) : 

oo oo oo 

Flix) = ^Pnx" + logx ^ + log^ X c„x", (101) 

n=0 n=0 n=0 



where the coefficients are given in eqs. fl94tl98p 

Let us stress that we have been able to obtain the complete analytical expression of 
the coefficients of this small-momentum expansion because of the knowledge of the MI at 
small momentum transferred. As we are going to show in the next sections, the absence of 
a similar knowledge in other expansion points will limit us to a numerical estimate of the 
corresponding coefficients. 

As expected from the threshold structure, the solution given in eq. fllOip is real for 
X > (space-like region) and has a non-vanishing imaginary part for x < (time-like 
region). The latter can be determined by simply giving the prescription for the logx factor. 
Since s ^ s + ie with e = -|-0, eq. (13T|) gives for x < 0: 

x-^-|x|-ie. (102) 

Consequently, we have: 

logx — >• log |x| — in ] 

log^ X — i> log^ |x| — vr^ — 2iTc log |x| . (103) 
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oo 



Therefore Fi{x) becomes complex for x < with: 

oo oo 

ReFi(a;) = ^p^x" + log |x| ^ g„x" + (log^ |x| - tt^) ^ c^x" ; (104) 

n=0 n=0 n=0 

(oo oo \ 

g„ a;" + 2 log |x| ^ c„ j . (105) 

n=0 n=0 / 

In the following sections, we will "move" on the real axis constructing the solution from 
X = to a: = oo; then with an analytic continuation we will switch to x = — oo and 
finally move back to x = through negative values of x. 

7.1.3 Improved Expansion — Bernoulli Variable 

A series for Fi with better convergence properties than the previous one in x can be con- 
structed by expanding in the Bernoulli variable [20] 



, / 1 + x\ 9 63 2 171 o 4095 . 32769 , / , , 

t = bg 8 = -X x^ + x^ x^ + x^ + O (x^) 106 

^\ 8 - xj 8 128 512 16384 163840 ^ ^ ^ ^ ^ 

= 1.125X - 0.492187x2 + 0.333984x^ - 0.249938x4 + 0.200006x^ + O (x^) . 
Eq. (11061) reahzes a one-to-one map between x G (—1, 8) and t G (— oo, oo). Furthermore, 

t = ^ X = 0, (107) 

while 

t ^ — oo for X — s> — 1"*" 

t +00 for X ^ + 8". (108) 

t diverges logarithmically when x approaches the singularities of the differential equation 
closest to the origin. Since this variable "follows" the singularities of the differential equa- 
tion, we expect a faster convergence with the order of truncation in t rather than in x. The 
inverse of eq. (I106p reads: 

= f:c.r. (109) 



X 



e^ + 8 

n=l 



The first few terms are 

28 o 44 , 35 . 1733 

9 

0.888888t + 0.345679^2 + 0.0603566 - 0.00533455 - 0.00586970 + O {t^) 



X = -t+—t^ + t^ + Oit^) (110) 

81 729 6561 295245 ^ ' ^ ' 



The radius of convergence of the above series is determined by looking at the singularities 
of X = x{t) in the complex t-plane, located in 

tk = 31og2 + i(2fc + l)7r, (111) 
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where k is an integer. The closest singularities to the origin, for k = 0,-1, give 



ro = y7r2 + 91og^2 ^ 3.76745. (112) 

We now substitute the r.h.s. of eq. fllOOp in the series expansion for Fi{x) obtained previ- 
ously and finally expand in powers of t, to 

oo oo oo 

F,{t) = F,{x{t)) = Y.^nt'' + logx + log^x J]7nt", (113) 

n=0 n=0 n=0 

where the coefficients /3„ and 7„ are determined from a^, hn and Cn respectively. The 
first few orders read: 

«o = -^ + ^; ^° = -]^' ^° = ^' ^114) 

"^= 648 - 72' ^^=108' ^^^ = -72' ^^'^^ 

3503 41 2 ^ _ 169 _ 41 

~ 2624400 ^ 11664 ' ~ ~ 87480 ' ~ 11664 ' ^ ^ 

We have computed -Fi(t) as a function of t (eq. flll3p ) by substituting the series fll09p into 
the series f llOip . The problem now is that of computing the radius of convergence po of the 
series expansion for Fi{t). The "internal" series for x = x{t) converges inside the circle 
in the t-plane of radius ro given in eq. (11120 : the "external" series for F{x) converges in 
the unitary circle |x| < 1 in the x-plane (eq. (!57|) ). As well known, power series always 
converge inside circles. The point is that a circle in the t-plane in general is not mapped by 
the function x = x{t) into a circle in the x-plane. Our problem has the following geometrical 
formulation: one has to find the largest circle in the t-plane satisfying |t| < ro, which is 
mapped inside the circle of unitary radius in x-plane by the transformation x = x{t) [21\ . 
On the unitary circle x = exp icp and one has to look for a minimum over ip of 



t 



1 + expiip 

log 8^ 

8 — exp lip 



16 

> log— < ro, (117) 



7 



where the minimum is obtained for ip = 0. We then findll: 
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po = logy ^ 0.826679. (120) 



^ Let us note that we did not express the logarithmic pre-factors log a; and log x as series in t, because 
there was no practical advantage for doing that. 

^ We have shown that the convergence radius of the t-series of Fi is not smaller than po, but it can 
actually be larger, depending on possible elimination of singularities. That can be illustrated with the 
following (rather trivial) example. Let us consider the expansion around x = of a differential equation 
having a solution of the form 

, , , , 8(1 + a;) , , 

0(x)=log^— (118) 

If we go to the Bernoulli variable t defined in eq. (|106p . we obtain 

4>it) = (j){x{t)) = (119) 
which can be analytically continued to all the i-plane, implying an infinite radius of convergence. 
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In general, for a series expansion centered in xq with nearest singularities a and b satisfying 



a < xo < b, (121) 

the Bernoulli variable t is defined as: 

f b — Xq X — a\ . . 

i ^ log ^ . 122 

\xo — a b — X J 

Note that 

t = <^ X = Xo, (123) 

while 

i — >■ — oo for X ^ a"*" ; 

t ^ +00 for a; ^ 6". (124) 
The following particular case are relevant in the following: 

• Xo — > oo: 

X — b J ' 

• 6 — > oo 

t ^ - \og( 

Xo - a 



^ = log(^); (125) 



t = -log(^^l; (126) 



• a — > — oo: 

t^-log( ^ ) . (127) 



XQ-b 



7.2 Expansion Around a: = 8 



In this section we consider the expansion of Fi{x) around the (space-like) point x = 8, 
which is a regular singular point of the differential equation not corresponding to any 
threshold or pseudothreshold of the diagram. The singular point of Fi (x) closest to x = 8 
is located in x = 0. We then expect the expansion around x = 8 to have a radius of 
convergence 

Rs = 8, (128) 
so that the series converges in the complex x-plane for 

|x -8| < 8. (129) 

For real x, that means: 

< X < 16. (130) 
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7.2.1 Homogeneous Equation 

By solving the indicial equation as in the previous section, we obtain a double zero in zero, 
so that the homogeneous equation has a solution of the form: 

oo oo 

fI^\x) = J]a„(a;-8)" + log(x-8)^6„(x-8)". (131) 

n=0 n=0 

The coefficients are, of course, different from those ones of the previous section — we use 
the same notation to avoid introducing too many symbols. The first few coefficients read: 

5 5 5 

ai^ -^ao - — bo; " ~ 24 ' ^^^^^ 

59 187 , , 59 , , , 

"^= 1728 "°+ 10368^°' ^^= 1728^°' ^'^'^ 

635 7561 , I, _ 635 , 

124416 "° " 2239488 ' 124416^° ' 

2171 59447 ^ ^ 2171 ^ 

04 = an H On ; Oa — On ■ 135) 

2985984 ° 107495424 ° ' ^ 2985984 ° ^ ^ 

We have again a triangular structure of the coefficients, as in the previous case. 

7.2.2 Complete Equation 

The expansion of the known term around a; = 8 is of the following form: 

oo 

C{x)= J^Qnix-Sr, (136) 



n=-l 



where the first three coefficients g„ are: 

= ^— 7r2 + — M0V2-— Mi + — M2V21og2 + — log^2; (137) 
^ 2304 768 256 192 ^ 256 ^ ' ^ ^ 

13 2 29 /- 13 29 ^ 1 , . 

TT^ MoV2H Ml M2V21og2H log 2 

82944 55296 ° 9216 13824 ^ ^ 768 ^ 

13 1 

log^ 2 log 3 : (138) 

9216 ^ 1024 ^ ' ^ ^ 

451 2 4637 /- 451 4637 /^r, ^ 

TT^ H MoV2 Ml H M2V2 log 2 

11943936 31850496 1327104 7962624 ^ 

551 , 451 2 61 , , 

log2 + log^2 + log3. (139) 

884736 ^ 1327104 ^ 147456 ^ ^ ^ 



18 



We have defined the following transcendental constants: 
Mo . r '^f^' + y^dy 

Jo Vvi^ + y) 

= — +41og21og(y2-l) + 21og'(V2-l)+4Li2 i(V2-l) +4Li2 -i{V2 - 1) 
b L J L J 

0.425435, (140) 

7^2 9 / 1\ 

Ml = — + -log' 2 + Lis f« 3.68568; (141) 



M2 = log(l + V2)^ 0.881374 . 

The general solution of the inhomogeneous 

00 



(142) 

> equation reads: 

Fi{x) = Y,^n{x - 8)" + log(a; - 8) J^Pni^ - 8)", (143) 

n=0 n=0 

with the first three terms given by: 

1111 1 
= -2304^^ - 768^°^+ 256^^ -192^^^^°S2--W 

-|ro-Ap„; (144) 

5 . ^ 

Pi = "24^°; (1^^) 

19 79 ^ 19 79 ^ 1 

r2 = tt' H MoV2 Mi H M2V2 log 2 log 2 

165888 221184 ° 18432 ^55296 ^ 3072 ^ 

19 , 2 1 , 59 187 , , , 

H log 2 H log 3 H ro H Po ; (146) 

18432 ^ 4096 ^ 1728 10368^ ' ^ ' 

59 

P. = YrasPo. (147) 

A particular solution can be obtained by setting po = 0, which makes all the coefficients 
Pn vanishing: that way the series multiplied by the logarithm disappears from Fi. One can 
also set ro = 0. The first three coefficients are then given by: 

ro = 0; (148) 
Po = 0; (149) 

ri = ^tt' - — M0V2 + — Ml - — M2V21og2- — log'2; (150) 

2304 768 " 256 192 ^ 256 ^ ' ^ ' 

Pi = 0; (151) 

19 79 ^ 19 79 ^ 1 

r2 = TT^H MoV2 Ml H M2\/21og2 log2 

165888 221184 " 18432 55296 ^ 3072 ^ 

19 1 

+ log' 2 + log 3; (152) 

18432 ^ 4096 ^ ' ^ ' 

P2 = 0. (153) 



19 



The general solution of the differential equation is given by 



Fi{x) = Ff^(x) + (154) 

and depends on the arbitrary constants ao and 60 entering fI^\x). Since we know from 
general arguments that x = 8 is a regular point for the solution, we can impose the 
logarithmic series be absent by requiring 



This condition gives: 



where: 



60 = 0. (155) 



F,{x) = J2^n{x-8r, (156) 

n=0 



So 


= ao ; 




5 


Si 








59 


S2 


1728 



(157) 

ao + n; (158) 
ao + r2, (159) 

etc.. Then, the expansion of Fi{x) around x = 8 given in eq. (11561) does not determine Fi{x) 
uniquely, because it still contains the free parameter oq. By using the matching procedure 
described in sect. 17.2.41 we obtain the following numerical estimate for this coefficient: 

ao = 0.0321062814000779405116. (160) 



7.2.3 Improved Expansion 

In order to improve the convergence of the series so far obtained, we move from the series 
in X to the one in the Bernoulli variable 

t = log|, (161) 

with the inverse: 

X = 8e' = 8 + 8t + At^ + ^t^ + + + ^t"^ + 0{f) . (162) 

3 3 15 90 

Unlike the previous case (x = 0), the above series has an infinite radius of convergence: 

rg = 00 . (163) 

Substituting the r.h.s. of the above equation in the series expansion for Fi{x) obtained 
previously and, finally, expanding in powers of t, we have: 

00 

Fi(t) = F,{x{t)) = (164) 

n=0 

20 



where the first three coefficients are: 



a2 



ai = 8 So ; (165) 

— ao + 7r^ + M0V2 Mi + M2 72 log 2 log 2 

54 " 5184 3456 576 864 ^ 48 ^ 

29 1 

+ ^log'2 + -log3; (166) 
343 203 o 929 , , /- 203 , , 929 



«3 = -T^cto - 



2 M0V2 + Ml M2y21og2 



486 46656 62208 5184 15552 

7 , „ 203 . 2^ 7 



log 2 log^2 log 3. (167) 

192 ^ 5184 288 ^ ^ 

A computation of the convergence radius pg of the series in t analogous to the one of the 
previous section gives: 

Ps = \og2 ^ 0.693147. (168) 



7.2.4 The Matching Condition 

oo can be computed by imposing that the series in eq. fllOip . which is completely determined, 
and that one in eq. f ll56p assume the same value in a given point in the intersection of the 
respective domains of convergence. One has to take a point lying in the interval 

- 1 < X < 1, (169) 
where the series of eq. fllOip converges, as well as in the interval 

< X < 16, (170) 

where the series in eq. fll56p converges. One has therefore to choose a point in the interval 

< X < 1, (171) 

such as for example x = 1/2. If we deal with infinite series, this procedure exactly de- 
termines the coefficient oq. As we have shown above, however, we can only determine an 
arbitrary, but finite, number of coefficients of both serieses and the matching has to be 
made in an approximate numerical way by using truncated series. The number of terms 
of the two serieses that must be computed depend on the required precision on oq. Our 
goal is to give Fi{x) with a relative precision of better then 10~^^ (double precision) using 
a relatively small number of terms in the serieses (around 30). The problem is that any 
point in the interval (11710 is close to the boundary of the convergence domain for the series 
centered around x = 8, where convergence is slow, implying that many terms are needed 
for high accuracy. In other words, a direct matching between the series in eq. (11 Oil) and that 
one in eq. fll56p is not the good algorithm. As we have shown before, a first improvement 
is obtained by re-writing the series expansions in terms of the relevant Bernoulli variables. 
But it is convenient to use additional series expansions centered around auxiliary regular 
points in the range < x < 8, such as for example x = 3. This topic will be discussed in 
more detail in sect. 17.51 
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7.3 Expansion Around a: = oo — Large Momentum Expansion 

In this section we consider the expansion of Fi{x) around x = oo, which generates a large- 
momentum expansion. Since the closest singularity to x = cx) is at a; = 8, we expect the 
expansion around infinity to be convergent outside the circle of radius 8, i.e. for 

|x| > 8. (172) 

The expansion at infinity is studied systematically by changing variable to 

(173) 



and considering the limit 

7.3.1 Homogeneous Equation 

Let us write the solution as usual as: 





1 


y = 






X 


y 


0. 



(174) 



F['\y) = S{y)Y,Any''. (175) 

n=0 

The pre-factor is assumed to have power-like form, 

S{y) = /, (176) 
and to be a solution of the homogeneous equation in the limit fll74p : 

S"iy)--S'iy) + ^Siy) = 0. (177) 

y y 

The indicial equation is (/? — 2)^ = 0, with a double zero in (5 = 2 and therefore the solution 
is, in the original x variable, of the form: 

oo oo J 

n=2 n=2 

Because of the presence of the logarithmic term in front of the second series, the infinity 
is a singular point — more exactly, a branch point of infinite order. Furthermore, since 
positive powers of x do not appear, x = oo is a regular point for the serieses above. Let 
us write a bunch of coefficients in terms of the lowest-order ones 02 and 62: 

03 = 2 02 - 3 62 ; 

33 

04 = 10 02 — ^ ^2 ; 

as = 56 a2 — 100 62 ; 

2561 

ttQ = 346 02 — 



^3 


= 262; 


(179) 




= 1062; 


(180) 


h 


= 56 62 ; 


(181) 




= 346 62 ; 


(182) 


h 


= 2252 62 • 


(183) 



4 ^' 

42631 , 
a-7 = 2252 tto bo ; 

Let us notice the power-like growth of the coefficients, which should scale asymptotically 
as ~ 8" (67/66 is already ~ 6.5). 
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7.3.2 Complete Equation 

The non-trivial part of the expansion of the known term C{x) around x = oo is related 
to the expansion of the GHPL's around this point, which is discussed in appendix |Bl The 
expansion is of the form: 



n=4 n=4 n=4 

where the lowest-order coefficients read: 

TT^ 7 

1 ISvr^ ,7 43 , , 

^^=2 + ^' ^^ = -4' "^^ = -16' ^^^^^ 

1 1257r2 , 131 331 

8 48 ^8 ^16 ^ ^ 

40 3577r2 , 3437 2569 , , 

= \ : I7 = : m-T = ; (188) 

^ 9 16 ' ^ 24 ' ^ 16 ' ^ ' 

18323 88277r2 , 56953 20301 , , 

kg = 1 ; <8 = ; = • (189) 

^ 288 48 ' 48 ' 16 ^ ' 

Since the differential equation involves a second derivative and the known term has double- 
logarithmic terms, the solution must contain up to four powers of the logarithm: 

00 00 00 00 + 

(-) = E ^ + i°g - E + - E + 'os' - E ^ + log' - E ^ • ( 190) 

•ij Ju Jb Ju Ju 

n=2 n=2 n=2 n=2 n=2 

Substituting the above form for the solution into the equation, we obtain the following 
relations among the coefficients: 

= 4; (191) 

S2 = 0; (192) 

= ^; (193) 
5 27 

P3 = 2p2-3g2--vr2 + y ; (194) 
vr^ 9 

gs = 2q2 + — + -; (195) 
16 8 

rs = (196) 
^ 16 48 ^ ' 

^3 = (197) 
ts - ^. (198) 
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The values of p2 and q2 are left undetermined by the differential equation and a particular 
solution can be found by imposing for instance 



P2 = 0; q2 = 0, (199) 
which we assume from now on. The general solution is given by: 

= Fi°\x) + Fi{x) (200) 

oo_ oo_ OO OO °° J. 

— + logx> — + log x> — +log x> — + log a;> — , 

rftlo ^^^^J rvi 1 1, ^^^^^ rt^ 1 1 ^^^^^ rp 1 1 ^^^^^ rp 1 1 

r, ^ O*^ O*^ O"^ O"^ 

n=2 n=2 n=2 n=2 n=2 

where we have defined: 

Pn = Pn + an] Qn = Qn + hn- (201) 

Let us note that the coefficients of the double, triple and fourth logarithm are uniquely 
determined because the associated homogeneous equation has no terms of this form: we 
therefore obtain exact analytic expressions for these coefficients. On the contrary, the 
coefficients of the terms with the single logarithm and of the terms without logarithm are 
determined up to a solution of the homogeneous equation, and then they do depend on the 
two arbitrary constants 

P2 and h- (202) 

The latter can be found by matching the values of the series in eq. (I200p with that one in 
eq. (1156P in a point in the range 

8 < X < 16. (203) 

By following the matching procedure outlined in the previous section and described in detail 
in sec. 17.51 we obtain 



P2 = -1.04850063265766512303; (204) 
g2 = +1.50257112894949285675. (205) 

By decomposing the two-loop integral into different infrared and ultraviolet regions, a large 
momentum expansion of Fi has been derived in ref. p^, which reads in our normalization 
containing an additional factor 1/16 : 



where ({3) = Yl'^=i ^/n^ = 1.20206 ■ ■ ■ . By comparing with the first-order terms of our 
complete solution, we find an analytical agreement in the coefficient of the fourth, triple 
and double logarithm|§. As far as the numerically determined coefficients are concerned, 
differences with the analytical expressions are at most O (10~^^), well within the expected 



^ We thank V. Smirnov for confirming a typo in the coefficient of the log^x term in [f 5', which must 
be divided by a factor 3. 
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numerical uncertainty. We can therefore use in place of our numerical estimates above the 
values: 

Our final expression for the large- momentum expansion of the amplitude Fi{x) of the 
crossed ladder diagram is that one in eq. (l2UUp . with the first coefficients given by: 

h - (208) 

m = +jC(3); (209) 

= (210) 

S2 = 0; (211) 

- f -|^'-lSo^'-T««'^ Pi^' 

9 TT^ 5 , , , , 

* = 8+16 + 2«<^)^ P"' 

13 TT^ ^ ^ 

= -l6-i8' <"^' 

= (216) 

<3 = ^. (217) 

Let us remark that eqs. (121311217!) provide analytic expressions of the coefficients of the 
leading power-suppressed terms in the large momentum expansion. 

7.3.3 Improved Expansion 

In order to accelerate the convergence of the series so far obtained, we move from the series 
in X to the one in the Bernoulli variable 

t = log^, (218) 
a; + 1 

with the inverse: 

1 1. 7 .2 11 3 , 35 . , 1733 . , 7217 . , 



^ — t 1' + + 1" + t'' + 0(f). (219) 

x 8 + e* 9 162 1458 52488 2361960 42515280 v ; v ; 

The convergence radius of the above series is equal to the one of the Bernoulli variable for 
the expansion around x = 0: 

7r2 + 91og2 2. (220) 
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Substituting the r.h.s. of the above equation in the series expansion for Fi{x) obtained 
previously and finally expanding in powers of t, we have: 

oo oo oo oo oo 

Fi(t) = 5^«„r + logx^/3X + log'a;5^7X + log'x5^M" + log'x^0„r, (221) 

71=2 n=2 71=2 71=2 71=2 

where the first three coefficients a„, 7„, and (pn are: 

«2 = ^aa; (222) 
P2 = ^h; (223) 

72 = (224) 
^2 = 0; (225) 
^ (226) 



15552 



5 9 11,5 , 

as = TT^ + 62 + a2 ; (227 

^ 34992 216 243 729 ^ ^ 

/Ss = —n^- — + —b2; (228) 

^ 11664 648 729 ^ ^ 

73 = —7r^ + ^^; (229) 

' 69984 11664' ^ ^ 



35 



139968 



343 



5038848 



where (see eq. (11991) ) 



The convergence radius of the expansion in t is: 



(231) 



95 2 2659 5 , 49 

04 = TT^ \ 62 H 02 ; (232) 

1259712 839808 2187 26244 ' ^ ^ 

5 29 49 

= " 104976'' 46656 + 26244''''' 

74 = -^^.^ + ^; (234) 

' 2519424 139968 ^ ^ 

35 

^4 = — (235) 

104976 ^ ^ 



(236) 



31 _ 5 

02 = P2 = tt"^: 62 = ^72 = +-;z(3). (237) 

2880 ' i ^ V ; v ; 



16 

= log— ^ 0.826679. (238) 
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7.3.4 Analytic Continuation to a; = — oo 

The expansion of the amphtude Fi{x) for large time-like momenta, namely for (cfr. eq. (11721) ) 

- oo < a; < -8, (239) 

can be found from the asymptotic expansion in the space-like region (x > 8) simply by 
analytic continuation. With the causal prescription 



X —\x\ — le , 

we have in addition to the replacements (11031) : 

log^x — » log'^ |x| — 37r^ log |x| — 3i7r log^ |x| + ivr'^ ; 

log'^x — > log'^ |x| — Gtt^ log^ |a;| + TT^ — 4^77 log^ |a;| + 4i7r^ log |x| 

Fi{x) then becomes complex for x < with: 



(240) 



(241) 



CXD ~ OO _ 



n=2 



n=2 



n=2 



(242) 



+ (log^ |x| - 3n' log kl) X] ^ + (W \x\ - log' \x\ + vr^) E ^ ' 



Im Fi (x) 



TT 



oo _ 



n=2 

oo 



n=2 



n=2 



n=2 

oo 



n=2 



(4W|x|-4vr^log|x|) J]^ 



n=2 



(243) 



7.4 Expansion Around x = —1 



In this section we consider the expansion around x = — 1, the only singularity of the 
differential equation in the time-like region. Since the nearest singularity to x = — 1 is the 
origin x = 0, we expect the convergence radius of the expansion around this point to be 



i.e. convergence in the circle 



On the real axis, that means: 



|x + 11 < 1. 



2 < X < 0. 



(244) 

(245) 
(246) 



7.4.1 Homogeneous Equation 

By solving the indicial equation, the solution turns out to be of the form: 

oo oo 

Fi'\x) = 5^a„(x + 1)" + log(x + 1) 5^6„(x + 1)". 



n=0 



n=0 



(247) 
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The differential equation allows one to express all the coefficients of the expansion above 
in terms of two of them, such as for example the lowest-order ones gq and 60: 

4 4 4 

01 = -ao + -60 ; h = -60 ; (248) 

41 62, ,41, , , 

02 = 27^0 + 81 ° ' ' " 27 ° ' 

400 2200, ,400, , , 

03 = an H On ; 03 = On . (250) 

243 2187 243 ^ ^ 

We have a triangular structure of the coefficients, as in the case of the expansion around 
X = 0. 



7.4.2 Complete Equation 

The expansion of the known term around x = — 1 is of the form: 



Cix) = J2 9nix + ir + logix+l)J2rnix + ir. (251) 

n=— 1 n=0 

The leading singularity for a; ^ — 1 is therefore a single pole and the subleading singularity 
is a logarithmic one related to the second series. The first few coefficients read: 

^" 4 216 1080 \72 540 / ' ^ ^ 

ro = ^; (254) 

2 131 o 1861 . /1627 1861 



- "3 + 972 



2 1861 /-^^ . /1627 1861 /-A 
47 

ri = ^; (256) 



Q2 



5177 1589 2 60178 /-^^ . /129527 120356 x 



4320 8748 273375 V 291600 273375 

749 

r2 = , (258) 

648' ^ ^ 

where K and i^' are two real transcendental constants defined as: 



K = f -j^M=^dy = - ^ fti - 1.97392; (259) 
Jo 



logy = _ ![! 
A/y(y + 4) ^ 5 



X' = log ( ilL^ j ~ 0.481212. (260) 
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Let us note that the coefficients above have an imaginary part related to the 2 massless 
particles threshold in s = 0. Unlike previous cases, the expansion point lies indeed in the 
time-like region. The particular solution of the inhomogeneous equation is of the same form 
as the homogeneous one: 

oo oo 

F,(x) = J2^r^(x + lr + \0g(x+l)J2tn(x+lr. (261) 

n=0 n=0 

For the coefficients, we find: 

si = ^so + ^to + q-i; (262) 

ti = ^to; (263) 

41 62 1 59 2 301 . / 7 301 /-\ 

S2 = — so + — to-- + + V5K + i7r{ VS ; (264) 

27 81 8 864 4320 V288 2160 / ' ^ ' 

1 41 

^2 = — + —to ; (265) 

16 27 ^ ^ 

400 2200 13 169 ^ 4381 r-^^ f 313 4381 , , 

= 243'^"+ 2187'°- M + 1944^ + 48600^"^ + ^^^ 
1 400 

^3 = - + to ■ (267) 

^ 8 243 ^ ^ 

etc., and we can choose Sq = to = 0. We assume therefore from now on the following 
expressions of the coefficients: 

So = 0; (268) 

to = 0; (269) 

^ "t^^K-C^i^K'- (270) 

24 24 12 ' ^ ' 

ti = 0; (271) 

1 59 2 301 . / 7 301 A 

S2 = — + 7r^ + V5i\: + «7r V5 : (272) 

8 864 4320 V288 2160 /' ^ ' 

^2 = ^; (273) 



S3 



13 169 2 4381 r-^^ . / 313 4381 /- ^ 



54 1944 48600 V 6480 24300 



^3 = ^. (275) 
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Finally, the general solution is given by: 



= Fi;°\x) + Fi{x) (276) 

oo oo 

= 5^5„(x+l)" + log(x + l)^t„(x + l)^ (277) 

n=0 n=0 

where 

Sn = cin + Sn and in = bn + 1„. (278) 
Eq. fl277p depends on the two arbitrary constants, 

Oo and bo ; (279) 

that can be determined in a numerical way as shown above: we equate the series expansions 
centered around x = and x = — 1, together with their derivatives, in one point where 
both series converge, such as for example x = — 1/2. That way we obtain: 

ao = -0.07572639563476980715 + ^0.3122156449500221544; (280) 
\bo\ < 10~^^ (281) 



naturally implying 



0. (282) 



Let us observe that for x < —1 the logarithm above acquires an imaginary part related 
to the 3-particle threshold at s = m^. Then, absorptive contributions related to both 
thresholds formally come out as imaginary parts of the coefficients of the serieses and, via 
analytic continuation, of the logarithmic prefactor. 

The vanishing of the coefficient bo of the log(l+a;) term in the expansion of Fi (eq. (I282p ) 
can be understood as follows. With a proper routing, Fi can be written as: 

V^kiV^k2 1111 1 

(g + kiY + ni? {q + + rn? {pi + ki)^ (p2 + {pi + ki - k2Y {p2 - + ^2)^ 

(283) 

By going to the threshold point = — m? and neglecting quadratic terms in the loop 
momenta ~ fcj, fci ■ /c2, one obtains in the soft limit: 

. / -i^ -1- -i^ -i^ (^«*) 

The integral above is convergent by power counting in the soft region. It is also easy to see 
that there cannot be any collinear singularity. To have coUinear singularities, one should 
have at least one "all massless" 3-point vertex, i.e. a 3-point vertex connecting massless 
propagators and/or external lines with light-cone momenta only, which is not the case. 
Since Fi has not any infrared singularity in the point = —m?, it is finite for x — >■ — 
and therefore cannot contain any log(l + x) term in the expansion, implying 60 = 0. 



Fi is ultraviolet finite, as shown previously. 
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7.4.3 Improved Expansion 

In order to improve the convergence of the series so far obtained, we move from the series 
in X to the one in the BernouUi variable 

t = - log(-x) , (285) 

with inverse: 

= _1 + + 1,3 _ J_,. ^ _L,= _ _L,e + o(f) . (286) 

The convergence radius of the series above is, as well known, infinite: 

r_i = oo . (287) 

Substituting the r.h.s. of the above equation in the series expansion for Fi{x) obtained 
previously and finally expanding in powers oi t, we have: 

oo oo 

F^(t) = F^(x{t)) = Yl + + 1) I] '^"^^ ' (288) 

n=0 n=0 

where the first three coefficients q;„ and /3„ are: 

tto = ao ; /5o = ; (289) 

a, = A + + IIL^K + t7r ( — - ^^/b] ; (3^ = - . (291) 

8 864 4320 V288 2160 /' ^^16 ^ ^ 

The convergence radius of the expansion in t is given by: 

p_i = log2 fti 0.693147. (292) 



7.5 Series Expansions Around Regular Points and Additional 
Matching Points 

In the previous sections we illustrated the method of the solution of the differential equation 
on the MI Fi{x), building scries around the singular points only. That is certainly necessary 
but not sufficient to determine Fi on the entire real axis. The expansions around all the 
singular points do not allow indeed to compute the master integral in the range 

- 8 < X < - 2. (293) 

To cover this region, it is necessary to perform expansions of Fi{x) around auxiliary (reg- 
ular) points Xi, expansions which are not multiplied by any singular function: 

oo 

Fi(x) = J2c^Hx-x,r. (294) 

n=0 
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It is clear that each new series expansion carries two unknown coefficients, which have to 
be determined by matching with known serieses, as we have made for the singular points. 
To cover the region (I293P , we have constructed two additional series expansions centered in 

X = -8 and x = -3. (295) 

Even though the space-like region < x < 8 is in principle covered by the expansions 
around the singular points, we have found it convenient to consider also series expansions 
centered in 

X = 1 and X = 3. (296) 
The matching procedure is the foUowinj^: 

1. we match the series centered in x = 1 with the series centered in x = 0, which is 
completely determined. That way we fix the two arbitrary coefficients entering the 
expansion in x = 1; 

2. we then match the series in x = 3 with the one in x = 1; 

3. we finally match the series in x = 8 with the series in x = 3. 

With the above method, we have been able to cover the range < x < 13 with a 
numerical precision of about 21 digits. In the range 8 < x < cxd, we need an additional 
series expansion in x = 16, to be matched with the series centered in x = 8. 

It is remarkable that the 9 expansion points above (4 singular + 5 regular ones) are 
sufficient to build a routine evaluating numerically the function Fi{x) for every value of the 
variable x in the real axis with a relative precision of about 10~^°. 

The internal consistency of the method can be checked as follows. The series expansion 
centered in x = is completely determined by the initial conditions that we have been able 
to find. Matching the free parameters of the series expansions centered in x = — 1 and 
X = 1 with the one centered in x = we obtain two other series without free parameters. 
Now we move along the real positive axis up to x = oo; 

x = — * x = l X = 3 X = 8 x = oo. (297) 

After that, we perform an analytic continuation in order to find the series expansion in 
X = — oo and, finally, we reach the series expansion in x = —1. What we find is a 
perfect agreement, within the required precision, between the free parameters that we fixed 
matching the series expansion in x = — 1 with the one in x = and those found performing 
the several matches along the real axis, back to x = — 1 0. Another (independent) check 
of our matching procedure is provided by the values of the coefficients in eqs. f l204ll205p . 
which are in agreement with the asymptotic expansion of -Fi(x) given in [15], as discussed 
before. 

The practical advantage of this method with respect to the general analytic continuation by (partially 
overlapping) circles is related to the use of the differential equation to generate both serieses. 

Since we deal in general with multivalued functions, one has to be careful to remain within the same 
sheet of the Riemann surface by using consistently the causal +ie prescriptions. 
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8 Expansions for the Master Integrals F2 and F3 

The second MI F2 is directly determined from the first one by means of the (algebraic) 
equation (l33l) with ^ = 1, which gives 

1 dF 

F2{x) = -x—^{x) + (298) 
/ ax 

in which Fi and F[ are assumed to be known. One just substitutes the power-series ex- 
pansions obtained so far for Fi{x) and obtains series representations for F2{x) in the same 
convergence domains. 

The third MI satisfies a first-order equation of the form 

^(x) + -Fs{x) + N{x) = 0, (299) 
where the known term can be written as: 

(300) 

It involves harmonic polylogarithms as well as Fi and F2, which are assumed to be knowij^. 
The associated homogeneous equation, 

^ + IfP{x)=0, (301) 
ax X 

has the general solution 

Ff (x) = -, (302) 

X 

where K is an integration constant. A particular solution of the inhomogeneous first-order 
equation is found in a straightforward manner with the general method of the variation of 
the constants: ^ 

F3(x) = - - f x'N{x')dx'. (303) 
X Jo 

Since Fi and F2 are not given in closed form in terms of known transcendental functions, 
but only in terms of (truncated) series expansions, one has also to expand all the GHPLs 
appearing in N{x). As in the case of the first MI Fi, the general solution is: 

Fsix) = - + F,{x). (304) 

X 

Analogously to the case of Fi, one can determine analytically all the coefficients of the 
expansion of F2 around x = 0, by fixing 

K = (305) 



It is interesting to note that formally Fi and F2 appear as sub-topologies in the evaluation of F3. That 
also implies that the general method of series expansions allows to compute topologies with subtopologies 
themselves given by series expansions. 
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in order to avoid infrared power divergencies. The matching procedure is similar to the 
one used for Fi, the main difference is that in this case there is only one undetermined 
coefficient for each series because the differential equation is of first order. Details can be 
found in 122|. 



9 Relation to the Equal-Mass Sunrise 

There is a conspicuous relation between the homogeneous differential equation for the first 
master integral Fi and the homogeneous equation for the equal mass sunrise in two dimen- 
sions, 



Siz) 



where: 



1 



D1D2D. 



(306) 



Di = kl + 

D2 = kj + m'^, 

D3 = {p-ki-k2y + m\ 



and 



(307) 
(308) 
(309) 

(310) 



with s = —p^ and the external momentum. The homogeneous equation (see eq. (4.3) 
of [13]) reads: 



^S(.;0) 



By writing 



1 1 1 

+ 



z z + 1 z + 9 J dz 



d 



-^(.;0) + 



3z 4(2+1) 12(z + 9) 



Fi{x) = -M{-x - 1) 



X 



whose inverse reads: 



M{z) = -{z + l)F^{-z-l), 



S{z) = 0. 
(311) 

(312) 
(313) 



we find that M{z) satisfies the same homogeneous equation of S{z). The "physical origin'' 
of such relation is unclear to us. 



10 Conclusions 

We have computed the three master integrals for the crossed ladder diagram with two 
equal-mass quanta exchanged by means of various power-series expansions centered around 
different points. The two-equal-mass exchanged non-planar topology is by far the most 
complicated one entering the two-loop form factor and a novel structure emerges, related 
to the occurrence of elliptic integrals rather than (generalized) harmonic polilogarithms. 
That was not expected a priori, because the diagram (see fig. 1) has only thresholds in 
s = and in s = m^, as well as a pseudothreshold in s = — 4m^; we may only generically 
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relate such structure to the non-planar topology. We argue that the occurrence of elliptic 
integrals is probably a general feature of massive multi-loop diagrams beyond some level of 
complexity. 

By studying the diagram close to the massless threshold s = we have obtained a small 
momentum expansion with coefficients analitycally determined. By combining our large- 
momentum expansion with that obtained in ref . [T5] , we have been able to give the leading 
power-suppressed corrections to the logarithmic contributions in a completely analytical 
form. 

We have found an a priori unexpected relation between the basic (first) master integral 
for two-mass crossed ladder in four space-time dimensions and the basic master integral for 
the equal mass sunrise in two dimensions: the homogeneous differential equation is exactly 
the same, while the inhomogeneous terms are unrelated (in the crossed ladder case, the 
latter ones being, of course, much more complicated). 

This work terminates the computation of the master integrals entering the two-loop 
electroweak form factor. 
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A Generalized Harmonic Poly logarithms 

The general theory of the harmonic polylogarithms has already been discussed in [TTl [TOl 
123] , to which we refer for details. Here we only give the specific information relevant to 
the computation of the crossed ladder. Let us first give the expressions of the harmonic 
polylogarithms (HPL's) entering our computation in terms of the standard ones: 




(314) 
(315) 

(316) 
(317) 
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Let us now show that the Generahzed Harmonic Polylogarithm (GHPlg 

iJ(r,0;x)= r —^M==dy, (320) 
Jo yjy{^ - y) 

is real for < a; < 4, purely imaginary for a; > 4, while it is complex for x < 0. In fact, for 
< X < 4 the integrand is real in all the integration domain and therefore also H is real. 
For X > 4, since x ^ x — ie with e = +0 and, therefore, y/ A — x + ie i\/x — A, one can 
split the integral as: 

H{r, 0; x) = f J^^^ dy - t f ^-^^ dy. (321) 
io Vy(4 - y) Ji y/y{y - 4) 

The first integral on the r.h.s. vanishes, 

logy 



\/y(4 - y) 



dy = 0, (322) 



implying that H is purely imaginary in this region. For x < the integrand is complex and 
so is the H. As a consequence, the term entering C(x), the inhomogeneous (and known) 
term in the differential equation for Fi given in eq. (l5l|) . 

H(r, 0: x) , , 

^ ' ' ^ (323) 



X 



is real for x > 0, as it should, while it is complex for x < 0. 

Let us now show that H{r, r, 0; x) is real for x > while it is complex for x < 0. Let us 
write: ^ 

H{r,r,0;x)= I ] H{r,0;y)dy. (324) 

Jo Vz/(4 - y) 

For < X < 4, the integrand is real in all the integration domain, while for x > 4 one can 
write: 

H{r,r,0;x)= [ j H{r,0; y) dy - t [ H{r,0;y)dy. (325) 

Jo ^/y{4: - y) A ^y{y - 4) 

Both integrals on the r.h.s. are real. 

H{r, 0; x) can be expressed in terms of ordinary harmonic polylogarithms of a non linear 
function of x, by changing the integration variable from y to 



« = I - ^ + (I - = f - 2 - J - ? - ? + « I • (™' 
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The weight "r" labels the integration over the function 



-''('^'*) = ^77^' (318) 

such that, for instance: 

, ^ r dt , , 

H(r;x)^ , (319) 

^ ' Jo Vti^ 
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where the square root is the arithmetical one for ?/ > 4, and on the last member we have 
made the expansion for y 3> 1. The inverse reads: 



it + l) 



(327) 



Note that y = 4 <-> i = 1. We obtain: 



H{r, 0;x) = -2iH 



X 



+ iH 



X 



In a similar way, we obtain for the only GHPL of weight 3 the representation: 



H{r,0,0;x) = -H 

+ 6C(3), 



0,0,0;--l + ^/a: 



(^-l) +4H 0,-l,0;|-l + Jxg-l) 



6 ■ 

(328) 



+ 
(329) 



where iJ(0, — 1,0;1) = —3/2^(3). In a completely analogous manner one can analyze 
H{r, r, 0; x), the GHPL entering the know term of the differential equation for F3. 



B Expansion of Generalized Harmonic Poly logarithms 

The points a; = 0, 4 and 00 are singular points for H{r,0; x). That is easily seen just by 
looking at the integral definition. The Taylor expansion of H{r, 0; x) centered around a 
regular point xq 7^ 0, 4, 00 is easily obtained by remembering that 

^H{r,0;x) = f{r,0;x), (330) 

where 

/(r,0;x) = ^^f=. (331) 

a/x(4 — X) 

By performing the Taylor expansion of /(r, 0; x) and integrating on both sides, one obtains: 

00 ^ 

H{r, 0; x) = H{r, 0; 0:0) + Yl Tu ^^"^ ' ^0)"^' " (^32) 

Let us note that one transcendental constant only is involved in this expansion, namely 
H{r,0; Xq). Let us now consider the expansions around singular points. It is clear that 
H{r, 0; x) does not possess an ordinary Taylor expansion around x — — the origin is a 
branch point of infinite order — where derivatives are not defined. This GHPL however 
has an expansion around zero involving semi-integer powers of x multiplied with up to one 
power of logx, 

00 oc 

H{r, 0; x) = x^/^ ^ + x^/^ logx ^ 1^ x", (333) 

n=0 n=0 
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where fc„ and In are coefficients and the serieses on the r.h.s. have a convergence radius 
R = 4 — i.e. up to the closest singularity. The representation above can be obtained by 
expanding around y = the regular part of the integrand, i.e. the factor l/\/4 — y, so that 

^ logyy-'/' + 4 log 1/1/'^' + 7^ logyy'/' + ^og yy"/^ + O {y'/^) 



^y{4: - y) 2 16 ""^^ 256 '"^'^ 2048 

(334) 

and using the result: 



/ 



logyy'-dy = —— logyy^^' - -——y^ + K (335) 
a + 1 (a + 1)^ 



The expansion around x = 4 is obtained in an analogous way, i.e. by expanding the regular 
part of the integrand around x = 4, log y/^/y, and integrating term by term. An expression 
of the following form is obtained: 

oo 

H{r, 0; x) = {x- (x - 4)", (336) 

n=0 

where 

go = 2 log 2; (337) 
1°S2 

91 = (338) 

92 = — — log 2. (339) 

^ 80 320 ^ ^ 

X = 4 is therefore a branch point of second order. Let us note therefore that the expression 
has no singularities for x — 4. 



By using the representation obtained in the previous appendix, we can directly obtain the 
expansion of the GHPL around infinity, which turns out to be of the form: 

oo oo ^ 

H{r, 0;x) = - i ^ - ^ ^ogx ^ - i - log^ x, (340) 

n=0 n=0 

where the two serieses on the r.h.s. are convergent for |x| > 4 and the lowest-order coeffi- 
cients read: 

0; (341) 

- 2; (342) 

- 3; (343) 

20 , , 

--. (344) 



ao = 


T' 


bo 


Oi = 


-2; 


hi 


a2 = 


3 

~2' 
20 


b2 


as = 




h 
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